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Based on an exact functional form derived for the three-point vertex function F, we propose a 
self-consistent calculation scheme for the electron self-energy with V always satisfying the Ward 
identity. This scheme is basically equivalent to the one proposed in 2001, but it is improved 
in the aspects of computational costs and its applicability range; it can treat a low-density 
electron system with a dielectric catastrophe. If it is applied to semiconductors and insulators, 
we find that the obtained quasiparticle dispersion is virtually the same as that in the one-shot 
GW approximation (or GoWoA), indicating that the GoWoA actually takes proper account of 
both vertex and high-order self-energy corrections in a mutually cancelling manner. 
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The electron self-energy £ is a fundamental quantity 
to control the quasiparticle properties in a many-electron 
system. Its accurate determination from first principles 
is recognized as a matter of central importance in many 
fields of condensed matter physics. In 1965, Hedin pro- 
vided a nonperturbative self-consistent approach to £ in 
a closed set of equations, relating £ with the one-electron 
Green's function G, the dynamic screened interaction W , 
the polarization function II, and the vertex function r. 1 ) 
This formally exact formulation, however, allows of no 
practical implementation in its original form, because we 
cannot calculate the electron-hole irreducible interaction 
/, a key quantity in the Bethe-Salpeter equation to deter- 
mine r, through its original definition using a functional 
derivative, / = <5£/<5G. Thus we are compelled to adopt 
some approximate treatments such as the GW approxi- 
mation (GWA) in which T is taken as unity. 

For more than two decades, successful calculations 
have been done for molecules, clusters, semiconductors, 
and insulators in the one-shot GWA (or G WoA), 2 ~ 7 ) 
but this is usually regarded as a too primitive approxi- 
mation, mostly because it is, in general, not a conserving 
approximation in the sense of Baym and Kadanoff. 8,9 ) 
In contrast, the GWA is a conserving one, obeying the 
conservation laws related to the macroscopic quantities 
like the total electron number. Upon implementation of 
this fully self-consistent GWA, however, we are led to a 
puzzling conclusion that the experiment on quasiparticle 
properties in semiconductors and insulators is much bet- 
ter described in the G W A than in the GWA. 10 ' 11 ) A 
similar puzzle is also found in atoms and molecules. 12 ' 

In metals, on the other hand, neither the Go Wo A nor 
the GWA works very well, 13 ) requiring us to include T in 
some way in treating systems possessing gapless excita- 
tions. Some schemes have already been proposed for this 
purpose, 14 ' 15 ) but they do not satisfy the Ward iden- 
tity (WI), an exact relation between E and T due to 
gauge invariance representing the local electron-number 



conservation. 16 ) In 2001, based on general consideration 
on algorithms beyond the Baym-Kadanoff one, 17 ) one of 
the authors (YT) proposed a scheme incorporating T in 
the GWA with automatically fulfilling the WI. 18 ) This 
GWT scheme (see, Fig. 1(a)) succeeded in obtaining the 
correct quasiparticle behavior in simple metals, but it 
encounters a serious difficulty in the low-density electron 
gas; convergent results for £ are not obtained, if its den- 
sity specified by the dimensionless parameter r s is larger 
than 5.25 where there appears the dielectric catastrophe 
associated with the divergence of the compressibility k 
at r s — 5.25 and concomitantly that of the static II in 
the long wave-length limit. 19 ' 20 ) Incidentally the GWA 
does not suffer from this difficulty, because n in it is not 
a physical one satisfying the compressibility sum rule. 

In this Letter, we provide a new exact functional form 
for r, based on which we modify the GWr into a scheme 
free from the difficulty originating from the dielectric 
catastrophe. In order to illustrate the power of the mod- 
ified scheme, which will be referred to as GWr^; (see, 
Fig. 1(b)), we show the results calculated for the elec- 
tron gas at r s = 8. If it is applied to semiconductors and 
insulators, we find that the quasiparticle dispersion self- 
consistently obtained in the GWTwi is essentially the 
same as that in the GoWoA, indicating that the GoW^A 
is superior to the GWA in the sense that for the sys- 
tems with gapful excitations, it actually takes proper 
account of the mutual cancellation between vertex and 
high-order self-energy corrections. This observation re- 
solves the above-mentioned long-standing puzzle on the 
GWA in comparison with the GoWoA. Here we empha- 
size that this cancellation is proved to be the case up to 
infinite order in an analytically rigorous way with clari- 
fying the assumptions needed in the proof, in sharp con- 
trast with the claims of a similar kind in the past; 21-23 ) 
they were inferred from the behavior of low-order terms 
in perturbation expansion for metals. 

Let us start with recapitulating the exact relations for 
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Fig. 1. Iteration loops to determine the self-energy (a) in the orig 
inal GWr and (b) in the improved scheme, respectively. 



systems with translation symmetry in which momentum 
p is a good quantum number. The Dyson equation relates 
G(p) with £(p) through G(p)- 1 = G (p) _1 -E(p) with p a 
combined notation of p, spin a, and fermion Matsubara 
frequency ioj n = «7rT(2n+l) at temperature T with an in- 
teger n. 24 ) The bare Green's function Gq(p) is written as 
Go(p) = (iu} n — e P ) _1 with e p the bare one-electron disper- 
sion. The Bethe-Salpeter equation determines T(p,p+q) 

by 

F(p,p+q) = l + ^2l(p,p + q;p',p' + q) 
p' 

xG(p')G{p' + q)T{p',p' + q), (1) 

where represents the sum ,X) p v With use 

of r(p,p+g), n(q r ) and E(p) are, respectively, given by 

H(g) = -^G(p)G(p+g)r(p,p+g), (2) 
p 

S W ^ -^G^+^^^r^p+g), (3) 

9 

with W(g) = V(q)/[1 + V(q)n(g)], where F(q) is the 
bare Coulomb interaction 47re 2 /q 2 . 

In Ref. 18, the concept of "the ratio function" was 
introduced to obtain an approximate functional form for 
T(p,p+q) satisfying the WI. By exploiting this concept, 
we have explored an exact functional form for T(p,p+q) 
and succeeded in obtaining the following form: 

T(p,p+q) = T^{p,p+q)T^(p,p+q), (4) 

where T^- a \p,p + q) and T^(p,p + q) are, respectively, 
defined as T^ a ' {p, p+q) = 1 — {I) PiP+q Il(q) and 

Gip+qy'-Gipy 1 



G (p+q) '-Gob) '-AE 



(5) 



p,p¥q 



Here an average of I, (I) p , p + q , and a difference in the 
self-energy, AE p pf9 , are, respectively, introduced by 



(I) P , 



P +q = -^2l(p,p+q;p',p'+q) 



xG(p')G(p' + q)T(p',p' + q)/Tl(q), (6) 



and 



/^ p ^ = Y J I{P,P+T,p',p l + q)[G{p' + q)-G{p% (7) 



as functionals of G and /. If I is exact, T^(p,p + q) 
is nothing but T(p,p + q) in Eq. (1), as can easily be 
seen from the very definition of (I) PtP + q , and AE p p+g is 
reduced to T,(p + q) — E(p), leading to T^(p,p + q) = 1. 
Thus Eq. (4) provides the same T(p,p+q) as that in the 
Hedin's exact theory. In reality, the exact / is not known 
and we have to employ some approximate /, in which an 
advantage of Eq. (4) over Eq. (1) becomes apparent; the 
former provides T(p,p+q) satisfying the WI irrespective 
of the choice of /, while the latter does not. 

Physically / takes care of exchange and correlation ef- 
fects in T(p,p+q) and it is well known that this physics 
can be captured by the local-field factor for the homo- 
geneous electron gas or by the Jastrow factor for in- 
homogeneous systems. In cither way, these effects arc 
well described in terms of a function depending only 
on the inter-electron distance, which justifies to assume 
that I(p,p + q;p' ,p' + q) depends only on q to write 
I (p, p+ q; p' , p' + q) = I{q). If this assumption is adopted 
in our exact framework, we obtain {I) p , p + q — I(q) and 
AE PiJ>+ ^ = 0. Then, by defining T WI (p,p+q) by 



Fwi{p,p+q) = — — i — - . ._■ 

Go(p+q) -G (p) 



(8) 



we obtain T(p,p+q) = [1 — I(q)H(q)]Twi(PiP+q), a result 
given in Ref. 18, leading to the GWr in Fig. 1(a). 

By substituting this result of T(p,p + q) into Eq. (2), 
we find that Tl(q) is written as 



n(g) = 



Rwi{q) 



l + I{q)U W l{q) 
with IV\vi{q)i defined by 

Kwi(q) = - ^2G(p)G(p+q)T WI (p,p+q). 



(9) 



(10) 



Then we can rewrite E(p) in Eq. (3) into 

= -^G(p+ g )M/( g )r W /(p,P+9), (ii) 

with W(q) = V(q)/{1 + [V(q)+I(q)}U WI (q)}. Combin- 
ing these results, we can construct the GWr^/ scheme 
shown in Fig. 1(b). This scheme is equivalent to the GWr 
in obtaining E(p), but it is free from the problem of the 
dielectric catastrophe, because it does not contain the 
calculation of H(q) inside the iteration loop. 

It also renders a great advantage to the reduction of 
computational costs to calculate 11(g) not directly but by 
Eq. (9) via Hwi(q), because Eq. (10) can be cast into a 
form convenient for numerical calculations as 

n(p+q)-n(p) 



e p+qf 



(12) 



where uj q is the boson Matsubara frequency and n(p) [= 
G(j))e lUn0 ] is the momentum distribution func- 
tion. Note that this expression very much resembles the 
one for the polarization function in the random-phase 
approximation (RPA) Ilo(g), which is given by 



Tl (q)^^Y, G o(P) G o(P + q) = J2 

P 



p+q ) 



-f(e P ) 



(13) 



t p+q" 
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where /(e) is the Fermi distribution function. 

Two comments are in order: (i) Since it is conserved 
on the microscopic level in our scheme, the electron num- 
ber is conserved on the macroscopic level as well. We 
can assure this conservation law by explicitly considering 
gauge invariance; because, as Baym discussed, 9 ) G trans- 
forms in accord with Go with the change of gauge, Twi 
is gauge-invariant, implying that the conserving property 
of S in the GWr^/ is the same as that without Ty^i, 
i.e., in the GWA. (ii) With use of Eq. (8) and the intro- 
duction of e p [=e p — J2 q W(q)/(iui q — Ep+q+Cp)], our scheme 
provides an integral equation to determine G(p) through 



( luJn -e p )G(p) = l + J2 — 



W(q )G(p+q) 



(14) 



On the assumption of I(q) — 0, this equation coin- 
cides with the one for obtaining the asymptotically ex- 
act G(p) in a neutral Fermi system such as the one- 
dimensional Tomonaga-Luttinger model 25 ) or higher- 
dimensional models with strong forward scatterings. 26 -* 
This coincidence clearly demonstrates the intrinsically 
nonperturbative nature of our framework. 

Basically I(q) is at our disposal; it can be determined 
either by perturbation expansion or by some nonper- 
turbative approach, but Eq. (9) suggests us to choose 
I(q) = — G+(q)V(q) with G+(q) the local-field factor. 
Note, however, that the meaning of G+(<z) here is differ- 
ent from the ordinary one that is defined with respect to 
IIo (q) instead of Tlwi (?) ■ Fortunately, we already know 
a good form for G+ (q) with taking account of this subtle 
difference, which is G s (q) in Ref. 27, satisfying the exact 
limit due to Niklasson 28 ) as |q| — > oo. 

With this choice of I(q), the GWTwi provides us the 
self-consistent £(p) in the electron gas for r s > 5.25, in 
spite of the existence of the dielectric catastrophe associ- 
ated with negative k. After analytic continuation (iuj n — > 
Ld+i0 + ) of E(p) to the retarded self-energy T, R (p,uj) with 
using the Pade approximant, we obtain the one-electron 
spectral function A(p,u>) [= — ImG fl (p, u)/ir}; an exam- 
ple is plotted in Fig. 2(a) at r s =8 and T = 0.001^ F with 
Ey the Fermi energy. The corresponding result for n(p) is 
given in Fig. 2(b), exhibiting a jump at the Fermi level, a 
typical Fermi-liquid property, though its deviation from 
no(p) [= 0(pf — |p|)] the step function with p-p the Fermi 
momentum is much larger than that at r s = 2, the typical 
density appropriate to many metals and semiconductors. 
We find an interesting result for the quasiparticle effec- 
tive mass to* at r s = 8; for |p| < 1.4pp, to* is larger 
than to c the free-electron mass, implying dominance of 
the correlation effect over the exchange one, while the 
opposite is the case for |p| > 1.4pp to give to* < m e . 
This crossover in to* never occurs for r s < 5 where to* 
is always smaller than m e . 18 ' 29 ^ 

In the crystalline case, each quantity involved in the 
GWTwi should be represented in the matrix form with 
respect to the reciprocal-lattice vectors {K}. For ex- 
ample, G(p) is a matrix composed of the elements 
{Gk.k' (p, *w n )} with p a wave vector in the first Bril- 
louin zone. For some quantities, we need to add the con- 
version factors transforming from the plane-wave basis 
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Fig. 2. (a) One-electron spectral function A(p,u) and (b) mo- 
mentum distribution function n(p) for the electron gas at r a = 8. 
For comparison, n(p) at r a = 2 is also shown. 



to the Bloch-function one in considering the matrix ele- 
ments; for example, IIo k,K' (<?) is given as 



n O K,K'0?) ; 



/(en'p+q)-/(enp) 



rm'per 



x(np|e- l ( ctfK )- r |n'p + q)(n'p+q|e i(ctfK ' ) - r '|np),( 1 5) 

where |?ip) is the Bloch function for nth band. 30,31 ) 

With this understanding, we have applied the GWTwi 
to semiconductors and insulators possessing a gap in elec- 
tronic excitation energies. Then, without detailed com- 
putations, the self-consistently determined quasiparticle 
energy E p in our scheme is found to be well approxi- 
mated by that in the GoWoA, as we explain in the fol- 
lowing. 

Let us assume that Tlwi(q) — ^oil) and I(q) = for 
the time being. Then, we may rewrite Eq. (11) as 

W (q) 



m = -E 



Goip+q^-Goip)- 1 



\{p)G{p)-\ (16) 



with W (q) = V(q)/[1 + V(q)U (q)j and \(p) a dimen- 
sionless function, defined by 



X(p) = X(p,iuj T1 



G(p + q)W (q) 



(17) 



The quasiparticle dispersion E p is determined by 
G R (p,E p )- 1 = 0, amounting to E p = e p + j: R (p,E p ), 
where we obtain the "on-shell" self-energy as 



£«(p,£ p ) = -]T' 



Wo(q) 



(18) 



by analytic continuation of T,(p) in Eq. (16). In deriving 
Eq. (18), we have paid due attention to the convergence 
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of X R (p,E p ) in gapful systems. In fact, provided that 
I(q) = 0, e p and the integral in the right-hand side in 
Eq. (14) are, respectively, reduced to E p and X(p), lead- 
ing to the behavior of G R (p, lo) for lo near E p as 

l + X R (p,E p ) 



G r ( P ,lj) 



E„ 



(19) 



For comparison, let us consider the self-energy in the 
G W a A, which is given by E (p) = -^2 q Go{p+q)W (q). 
By analytic continuation ioo n — >• e p +i0 + , we obtain 

W (q) l w , , 



S?(p,e p ) = 



E 



9 e p+q + e p 



coth 



e p+q 



2T 



tanh 



e p+q 
2T 



(20) 



Because the transition p+q — > p involved in Eq. (20) 
is relevant only for the interband transition, |ep+q — e p | 
is always larger than E g the energy gap. At low T, the 
chemical potential fj, lies at the center of the band gap, 
indicating that |ep+q| >E g /2. These two facts allow us to 
safely neglect the contribution from the second sum in 
Eq. (20), as long as T <C E g . Thus we may write E p the 
quasiparticle dispersion in the GqWqA as 



£ p =e p + ££(p,e p )=e F 



E 

g 



W (q) 



(21) 



leading us to conclude that E® = E p . Note, however, 



that the spectral weight 



(i-d^(p,u)/duy 



is different from l + \ R (p, E p ). 

In the literature, E® is sometimes evaluated as E® = 
e p + z p S^(p, e p ) and there is a controversy as to whether 
this z p should be included or not. As previously discussed 
in detail, 23 ' we consider it better not to include z p so 
that the vertex corrections beyond the RPA are properly 
included, together with higher-order self-energy terms in 
a mutually cancelling manner. In fact, our present result 
of Ep = E p without this factor z p indicates that this 
feature of mutual cancellation reaches far up to infinite 
order in semiconductors and insulators. 

Finally we comment on the two assumptions: (i) The 
difference between TLwi and Bo arises only from that be- 
tween n(p) and nn(p). In usual semiconductors and insu- 
lators, the valence-electron density is high; for example, 
r s = 2 for Si. Now n{p) in a metal at such r s does not 
deviate much from no(p) except for the states near the 
Fermi level, as shown in Fig. 2(b), but those states are ab- 
sent from the outset in these gapful systems. Thus n(p) 
is close to no(p), leading to U WI w B . (ii) Justification 
of I = has already been done by numerical studies in 
Ref. 3, in which I in our scheme is critically assessed in 
terms of K xc the density-derivative of the Kohn-Sham 
exchange-correlation potential. From an analytic point 
of view, it is enough to note that the basic processes to 
contribute to I are related to the interband electron-hole 
interactions, in which |q| for principal processes is of the 
order of |K|, making V(q) very small and G+(q) reach 
its asymptotic constant. Thus the effect of I is weak in 
semiconductors and insulators. 

In summary, we have proposed the GWTwi scheme 



for the fully self-consistent and conserving calculation of 
the electron self-energy. This can be applied not only to 
metals in a wide range of densities but also to semicon- 
ductors and insulators, in which the obtained quasiparti- 
cle dispersion is close to that in the GnWnA, explaining, 
from a fundamental viewpoint of many-body physics, the 
reason why the GqWqA better describes the experiment 
than the GWA in those gapful systems. We also realize 
that the role of I, representing short-range exchange and 
correlation effects, is very much different between gap- 
less and gapful systems; in the former, it can never be 
ignored to obtain the reliable quasiparticle behavior, but 
in the latter, it can be neglected, as long as the GnWoA 
well reproduces the experiment. In this respect we can 
suggest that the Go Wo A should be performed with judi- 
ciously choosing the basis functions to make I and the 
difference between Uwi and B as small as possible. 
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